


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1980 


Optimal control system design with prescribed 
eigenvalues via Cauer Second Form. 


Stanley, Edward J. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/19003 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


f (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


OPTIMAL CONTROL SYSTEM DESIGN 
WITH PRESCRIBED EIGENVALUES 
VIA CAUER SECOND FORM 


Edward J. Stanley Jr. 





NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





Peo 


OPTIMAL CONTROL SYSTEM DESIGN 
WITH PRESCRIBED EIGENVALUES 
VIA CAUER SECOND FORM 
by 


Edward J. Stanley, Jr. 


September 1980 





Thesis Advisor: M. J. Goldman 


Approved for public release; distribution unlimited 


1197051 











SECURITY CLASSIFICATION OF THIS PAGE (Whan Data Entered) 


REPORT DOCUMENTATION PAGE 


2. GOVT ACCESSION NO 





READ INSTRUCTIONS 
ihe a COMPLETING FORM 


$. TYPE OF REPORT &@ PERIODO COVERED 
Master's Thesis; 


september 1980 


6. PERFORMING ORG. REPORT NUMBER 







4. TITLE (and Subtitie) 

















Optimal Control System Design 
with Prescribed Eigenvalues 
via Cauer Second Form 








AU THOR(8) 


“Edward i stanley, Jr. 


- CONTRACT OR GRANT NUMBER ea) 













ENT, PROJECT, TaS« 
NIT NUMBERS 







9. PERFORMING ORGANIZATION NAME AND ADORESS 


Naval Postgraduate School 
Monterey, California 93940 














12. REPORT CATE 


september 1980 


13. MUMBER OF PAGES 


- MONITORING AGENCY NAME & ADORESE/(I! differant trem Centrailing Office) 18. SECURITY CLASS. (of thte report) 


1t CONTROLLING OF FIZSE NAME ANO ADDRESS 


Naval Postgraduate School 
Monterey, California 93940 











Unclassified 


$e. O€CL ASSIFICATION/ DOWNGRADING 
SCH EOULE 


16. OISTRIBUTION STATEMENT (of thie Report) 


Approved for public release; distribution unlimited 


17. Ost RiIBUTION STATEMENT (ot the enatract entered in Bleck 20, ti ditterent trem Repert) 


18. SUPPLEMENTARY NOTES 


18. KEY WORDS (Continue on reverese eide it neceeeary and tdentity by block number) 


Cauer Second Form; Continued Fraction Expansion; Continued 
Fraction Inversion; Presceibed Eigenvalues; Transformation 


from State to Phase Variable Form 





20. ABSTRACT (Continue on reveree etde if neceeeary and identify by bieck number) 
A method is developed in terms of the Cauer Second Form repre- 
sentation of continued fractions as a means of designing linear 
Single-input-output (SISO) control systems. Optimal closed loop 
solutions corresponding to a set of prescribed eigenvalues are 
obtained through minimization of a quadratic performance index. 
The partitioning method of the Cauer Second Form for system sim- 
plification is presented with a simplified inversion technique 
for the reduced order system. 








DD i ani es 1473 EOITION x: {NOV 6613 OBSOLETE 


(Page 1) S/N 0102-014- 6601 : 
] SECURITY CLASSIFICATION OF THIS PAGE (Phen Deta Entered) 





Approved for public release; distribution unlimited 


Optimal Control System Design 
With Prescribed Eigenvalues 
Via Cauer Second Form 


by 
Fanwcdeceummouanley, or. 


Captain, United States Marine Corps 
Bec ev atlaneva University, 1972 


Submitted in partial fulfillment of the 
requirements for the degree of 


Peete ot SCIENGE IN ELECTRICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
September 1980 





ABSTRACT 


A method is developed in terms of the Cauer Second Form 
representation of continued fractions as a means of designing 
linear single-input single-output (SISO) control systems. 
Optimal closed loop solutions corresponding to a set of 
prescribed eigenvalues are obtained through minimization of 
a quadratic performance index. The partitioning method of 
the Cauer Second Form for system simplification is presented 
with a simplified inversion technique for the reduced order 


system. 
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Peet LT RODeC rT LON 


The purpose of this research was to develop an algorithm 
for obtaining optimal closed loop solutions corresponding 
to a set of prescribed eigenvalues for single-input single- 
output (SISO) control systems. It was desired that the 
algorithm be adaptable to digital computer techniques and 
unrestricted by system order. 

The Cauer Second Form for system dynamics representation 
was chosen over other alternatives because of the regular 
pattern of the state and output matrices, and the method 
of linear system simplification. 

In Chapter II, several basic properties of both Cauer 
First and Second Forms are presented from the theory of 
continued fractions. A simple and efficient algorithm is 
also developed for inversion of the continued fraction in 
either form, independent of Routh's algorithm. 

In Chapter III, the method of linear system order 

reduction based on the Cauer Second Form is amplified. 
The emphasis on this area was primarily to elucidate the 
various methods previously employed for system simplifi- 
Gation. 

The original theoretical work of this thesis is pre- 


sented in Chapter IV. The objective was to obtain closed 





loop solutions corresponding to a prescribed set of eigen- 
values. While minimizing a certain cost function, which 
met desired system characteristics. It is shown, by examples, 
that the derived algorithm is equally capable of handling 
systems with multiple and/or complex, as well as, unique 
sets of real eigenvalues. 

The final chapter, Chapter V, presents a discussion of 


results and suggests areas for future study. 





ees wore nRi eo OF CAUERK FIRST AND SECOND FORMS 


fees oeD LOOP SYSTEM IN CAUER FIRST AND CAUVER SECOND 





BORMS 
Consider the closed loop transfer function given by: 

n-l ; 

a b.s- 

mcs) _ L500 
> , onl ay ¢2-1) 

Ste a.S 

1=0 


with block diagram as given in Figure 2.1. Equation (2-1) 
can be expanded into the Cauer Forms of continued fractions 
as follows. 
1. Cauer First Form 
a. Arrange the numerator and denominator poly- 


nomials in descending order. 


bee rertorm cOonelnued division. 

n-l n-2 
Y(s) _ b-18 55 b-28 es ot bjs + Do 3 
UCs. sh +a sl 4 4 ge ee 

n-1 n-2 IL O 
= 1 

h,s zt 1 
lat il 
2 i a (2-3) 
h,s - ds 
hy f 


ag 
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2. Cauer Second Form 


a. Invert the numerator and denominator and arrange 


the polynomials in ascending order. 


n-2 n-L ia 
wey to «CF ee ae v ea. 
Us) ~ n-2 n=-1 a 
BG z b,s Te Paths a MME Ore c: <-> tD__os 5 By 
b. Perform continued division. 
mes) _ 
Wis) ~ i 
AL + ; 1 
2 
=< was 
= ie (2-5) 
h. i 1 
hes 
or 
ib 
lal ag S 
A a OE) 
+ S 
Ao 
he - S 
ny a 


Block diagrams of both systems are shown in Figures 2.2 and 


2 oe 


eZ 
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PeeeeeHYSICAL INTERPRETATION OF DOMINANT TERMS 
RESULTING FROM CONTINUED FRACTION EXPANSION 


It is known that the most dominant terms in equations 


(2-3) and (2-5) are the first quotients, h,s and h,>. re- 


1 
spectively. A meaningful interpretation for these terms 

can be accomplished by applying the initial value and final 
value theorems.” Letting Y(s)/U(s) = Fs), by an asymptotic 


expansion approximation: 


a. For Cauer First Form 


Lim £(t) > lim sF(s) = = (2-7) 
il: 
to S->0o 
and 
inne Ce) tam sFCs) ~ h, + h,. (2-8) 
too S>+O 
meee FOr Cauer Second Form 
mm £(t) ~ Lim sF(s) * he oF hy (2-9) 
t7o So 
teen f(t) > Jim sF(s) ~ hi: (2-10) 
tc S>O 


+ : F 
fam £(t) must exist. 


t-¥o 


is 





Equations (2-7) and (2-10) are of considerable interest since 


they involve the dominant term, h The implication is that 


1° 
the Cauer First Form emphasizes the initial or transient 
response of the system; whereas, the Cauer Second Form em- 
phasizes the final or steady state response of the system. 
In general, the quotients lower in position in the continued 
fraction expansion have less influence on the performance of 
the system as a whole, ane has less influence than ns, where 
1<j}). Because many systems must meet a set of steady state 


conditions, the Cauer Second Form will be used for the 


prescribed eigenvalue problem. 


wee CONTINUED FRACTION INVERSION 

The theory of continued fractions was first associated 
with Routh's Algorithm Deven tos. ti} and 12). The 
following year Frank [3] extended and modified Wall's work 
to include complex coefficients. Both, however, applied 
feuen's algorithm only to continued fraction expansions, 
not to the problem of inversion. 

In 1969, Chen and Shieh L4] developed an algorithm 
method for converting a continued fraction into a rational 
fraction of two polynomials. Their method, which makes use 
of Routh's algorithm, is presented below. 

If the elements, h;, are known for any continued 
fraction, then the state and output equations can be 


written immediately from Figures 2.2 or 2.3. 


ls 





Z4 hoh 
Zn hph 
Z 2 hoh 
Zn hah 
= 
al 
a8 
Zz 
n 


ois + DY 


Th, h 


h 


h, Ch, th.) 


h, Ch, +h.) 


h, Ch, +h.) 


ae 


Hiei Sees 


6 i Tien 
he Ch, +h.) +» hy Ch, ths) 
he (h,+h,+h,)-.. 
he Ch, th,th,)..-h, (Ch +--. thy, 4 
(2-11) 
Oi) 
ei 
ao 
(22D) 
Z 
nN 
(2-14) 





The coefficients of the characteristic polynomial, 
|sI-H|, become the first row elements of the required Routh 
array. The next sequence of steps in determining the (j+2)th 


row of the Routh array is to successively let 


and 


h. = 0 (2-15) 


fuego l,3,5,...2n 1) 
and evaluating the remaining a oT) a) matrix, where 
feet) )/2, up to k=n-1, i.e., for n arbitrary and j=1; 


the 3rd row of the Routh array becomes (after h. and h, are 


i 2 
set equal to zero) the coefficients of 
| sI-H, | ) 
where the (n-1)x(n-1L) Hy Measkiet x os: 
hh, heh. GENS eT oS aN Ahh 
~ thoh. h.(Ch,th_) h,_(h.th-.) omron 
y 3 Ge ais aoa Z aS 
Mee H.C th.) ........ nen Gaeya ieee) 


This process repeats until the system state matrix is 


reduced to a single element, H yielding the (2n-1)th 


2n-1’ 


row in Routh's array. It is observed that each successive 


ik 





odd numbered row contains one less element than it's 
predecessor. By inserting leading zeros in the 3rd, 5th, 


Mmeeentrl!)th row, the matrix, P, is formed. 


3rd eee Paez P13 

5th 0 eZ P23 

mtn 0 0 Es CZ ) 
mene) th 0 0 0 





The matrix, es 1s the linear transformation matrix required 
to obtain a linear system in Cauer Second Form from phase 
variable (canonical) form. Continuing, the second row of 
the Routh array is obtained from the output matrix, L, 


and the above transformation: 


ce= LZ Ceontl mleadmtractioens) 
w= Px (linear transformation) 
y = Cx (output equation, phase variable form) 


Therefore, 
@e= i P . (2-18) 


C is an (lxn) vector whose elements are the seoond row of 
the Routh array. 
Consider the Routh array as an (ntl)x(2nt+t1l) matrix with 


typical element ee The quotients, he of the continued 


IE 





fraction expansion can be expressed as: 


aul 
—— Geilo 


From this relationship and knowledge of how the Routh 
array is generated, the remaining even numbered rows of the 
array can be found. The transfer function as a ratio of two 


polynomials is written as: 


noms 
ry 


oe 


T(s) = i (2-20) 


| 
ih + 
2: 


Os 
fA 
fj 
i 


Chen and Shieh [4] contend that this method is the 
easiest in attaining the inversion. The author disagrees 
and presents a simpler iterative method based on the in- 
version technique for the Generalized Cauer Form given by 
Goldman [5]. The method is equally suited to both Cauer 
First and Cauer Second Forms, requiring no prior knowledge 
of Routh's algorithm. Assuming all h.'s are known, and 


non-zero, in equation (2-3) or (2-5), let: 


mee = Nl. - (2-22) 
aL Zale 
epee? G61,2,....,nJ. 


20 





lee Lnversion of Cauer First Form 


iphialairzZeuuvoOucnrixil) vectors C and D. 


©) 
I 

= 
Q 


1 co ttt cw] (2-23) 


O 
is 
a 
O. 
A, 
A. 


D 9 dy dg cess dl (a= oin) 


to all zeros, except: 


Ch = bo (2-25) 
d-i =a, xc (2-26) 
dq, ool. (2-27) 


The following set of equations are first solved for i=l. 


os. = ; Bee ty des =< 
et bi * See eee (2-28) 

; . os ‘ oh ase ; é - 
dn (itl) +4 eed x n-it3 aoe (atl) +4? (2-29) 
where j£L0,1,2,...,iJ are substituted in ascending order, and 


(2-28) 1s solved before (2-29) for each value of j. Now, 
let 1=2 in equations (2-28) and (2-29) and repeat the same 
procedure. The index "i" is incremented until i=n-1l, and 
(2-28) and (2-29) are solved as before over the appropriate 
range of the index "j". The final vectors, C and D, contain 
elements which are the coefficients of the numerator and 


denominator polynomials, respectively, of the transfer 


21 





fmemeclLon (or driving point impedance function): 


a n-1 
ees 
ccm 
ee (2-30) 
n=-J 
las 
jso 7 
Eycemp le : 
2 
cc 
T(s) = “gs ti7ist360 (2-31) 


Spe es en O2S+./ 20 


By continued fraction division: 


10s2#171s+360 |s?+71s2+702s+7206 .1s 


aaa sc 


53.85-+666s+720 


53.95°+666st720 hos?+171Ss+360 og meio 


ee Se eta 
) n7.438st226.u2 


47 43854226 .42 [53.9s2+666s+7204 1.1362s 


53.98°+257. 258s 


408.742s+720 | 


408.742s+720 ; st ° mS 
47 .438s+82.79 
143.63 (2-32) 
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143.63 fsO8-7H2S+7 201 2.8458s 
4O8.742s 


20 


720 pus. 63 AIRS Sis ) 
143.63 


0 


Therefore, the transfer function in the form of equation 


(2-3) is: 
T(s) = 1 
eS et if 
ee 5 524+ 1 
Leiso2s. + alk 
poli dc Oreect ah 
2.0455 S °r — | 
21995 
(2-32) 
wath 
hy =e Ay =) lS 
Ao eo 5S he Se EUS TS 
he =e i 2 he = wooo. (2-33) 


Now, using equations (2-21) and (2-22); 


ae 





a= h, = .l Peat = ALO OS 


il as alt Z 
ay = h Sere bs OFZ De = hy Meee le (2-34) 
az = Ae = 2.8458 D, = Ae =— 2 7h99 5 
From equations (2-25) through (2-27), 
Ca = b, =e. 
d, = A,xC. = eee Gs. 905) = «560 C235) 
d, = 1 


Substituting i=l into equations (2-28) and (2-29) for j=0: 


Cy = D,xd,tc, ecoCooo). + 0 = .0653 

dy - a5xCotd, —meliCoZ Uoos) + 0 = .0742 
for j = 1 

C2 = Bocd,tc, eeeio Gd) +t 21995 = 4.3145 

d., ~ a,xC,+d, TierOne lot. «568 = «9353 


At this point, j=1, therefore increment index "i": 


mem 1=2,  )=0: 


1 By xd, tc, =eeoe oC. 0/42) + 0 = 20138 


Q 
i 


0 a,xc, td, e038 + 0 = 7 00136 


ay 
i 
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Cy = Db xd,tc, See o eos) 0653 = .2388 

ds = a, xcytd, see soo) f 70742 = .0981 
for j}=2: 

C4 = Bb xd.tc, =o) to 23'S = . 5000 

d, = a,xc,td, —enGz oy + 90353 = 2.9853 


Now, at the point where j=1, and i=n-1, the transfer 


mane tion is: 


2 
e ay fe 
1 (2-36) 
POOmses +t. UG8IS +.9853stl 


Multiplying numerator and denominator by 1/d. = 20m yLTelas: 


2 
+ 
T(s) = yes =A tisteee _ Casi 


Say lsat /O2st7 20 
2. Inversion of Cauer Second Form 
Initialize the two (ntlxl) vectors C and D, (2-23) 


and (2-24), to all zeros, except: 


= bE (2-38) 

do-1 aes (2-39) 

ae = a xc (2-40) 
nh Ti Ti 


295 





The following set of equations are first solved for i=l. 


4n-(i41)45 * An-ad * Cn-citaej * Sn-it1 (2-42) 
: sate ‘ : : 
where jEL0,1,2,....,;i1) are substituted in ascending order, 
and (2-41) is solved before (2-42) for each value of j. 
Next, find at eceond ing.) Ee! 
m= as XC. (2-43) 


Now, let i=2 in equations (2-41) and (2-42) and repeat the 
Same recursive procedure. The index "i" is incremented by 
one until i = n-1, and for each value of i, (2-41) and (2-42) 


are iteratively solved over the appropriate range of the 


ro Tt 
e 


index "3 The resulting elements of C and D are the 
coefficients of the numerator and denominator polynomials 


Sumehne transfer function: 


(2-30) 


Using the sample example, (2-31); 


10s*+171s+360 


5°4+7197+7028+720 


TCs) = 
+ 
fem j=l, c 


n-itjt+l = Q. 
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Place the numerator and denominator terms in ascending 


order, invert, and perform continued fraction division: 


360+171stl0s~- /720+702st71s ‘+s l 2 


720+3425+205° 


epee tee 
3 2 
360st5lsts [360+171s+10s L l/s 
360+51sts- 
120s+9s" 


120s+9s- [36054515 +50 3 


360s+27s- 
2s “+s° 


ets [120+9st 5/s (2-44) 


iP2oSs 
US 


24u+s 


= 

N 

= 
mo 





9) 


s [tl 4/s 
i 
0 
The transfer function (2-31), in the form of equation (2-6) 
aS 
L+ S 


3+ S 


Zi 





a 
i 


5s 
il 


ar 
Tt 


oS 


Z hy 
i he 
3 he 


S 
Bt S 
ats) 
uy 
= 5 
= 6 
= lit (2-46) 


For the inversion process, using equations (2-21), (2-22) 


and (2-46); 
a, = hy = 2? by 
an = h. = 3 Db, 
a, = A. = 6 b. 
and from 


© 
it 


= he =~) 
=h, = 5 
=h, = 4 Qa 


equations (2-38) through (2-40), 


b, = 4 
aE 
a, xX C, = 6(4) 


= 24. (2-48) 


Now, substituting i = 1 into equation Oe a sande 2-42 ero 


}=0: 


3(0) + 1=1 (2-49) 


Ne 





hor 


and 


for 


for 


OY 


and 


b, Xx d, + O) S SC) aibate 
a, x C, + d, sop 24 = 5) 
equation (2-43): 
an x C 3 = 3(120) = 360 
j=0: 
by = Co =l]1+g9= 10 
(2-50) 
a, X cp + ds = 2(0) t+1letl 
by Xx d,, + C3 = iC Si) os Ibe Sb 
: = + = 
a, x cy + d, PAA Sak aol 
by x d. ne = eles 60) eaten = 3.5.0 
a, xX Cy zy d. ——JGle yap e some = 7012 


equation (2-43): 


Been CG r= 2.5900) = 720. 


3 


Zo 





Therefore, 


LO 10 tc 360] 


Ce 
i 


De— {1 ol 702 ae. (2-51) 


and the transfer function realized from equation (2-30) is: 


Z 
re 
T(s) = ti tg fSteet CRE) 


Sea hs t/702s+ 720 
which is the same as (2-31). 

This completes the development of the continued fraction 
inversion algorithms from Cauer First and Second Forms. This 
iterative procedure is easily seen to be computationally 
much simpler than Chen and Shieh's method. First, it does 
Gemeeequire the need to find the H matrix, (2-11); and 
second, it does not necessitate finding the coefficients 
of n characteristic polynomials of diminishing order. This 
method is solely based on equation (2-6), enumerating the 
inversion from bottom to top. As by-product, the entire 
Routh array appears in the intermediate steps as can be 


seen from the Cauer Second Form example: 


eo, 64, «(d, Ome O2 oy wal) 

mc, Cy 360 171 #210 

ad, d, = 360 51 1 

oo 120 g 

ac 24 x 

7 . 

1 i (2-53) 
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where rows 2(n-i)1 and 2(n-1) are taken from the ith iteration, 
Mngt 2,..-.,n-l}]. "“"i=0'' implies the rows come from the 
ieGialization of C and D. The last row, the C7nety ths 2s 
always the single element one. 

If a comparison is made between equations (2-20) and 


(2~30), it is observed that: 


i na +1” 1€L1,2,...,nJ 


ae = 


5 1,x-j+1? Repay 2 tere ent 3 (2-54) 


where the c.'s and do's are taken from the (n-1)th intera- 
tion under the index "i". 

It is also observed that if the quotients, h.'s, 
resulting from expansion into Cauer Second (First) Form 
are used in the inversion algorithm presented for Cauer 
Farst (Second) Form, then the Cc. 's and d.'s in the (n-1)th 
iteration represent the transfer function coefficients in 
reverse order. This is shown using the preceding example. 


From equation (2-46); 


h, = 2 h, = 5 
h, = 1 nee 
me = 3 he = 4 


and equation (2-47); 


a, =h ee b, =h = 1 
ao = bh — ae bo = hy = 
a, = he = 6 Db. =h, = 4 





Now, using the inversion scheme for Cauer First Form, from 


(2-25) 


through (2-27); 


xO. = 6(4) = 24 (2-55) 


5 & Gite e= 54 24) 4.0 =. 120 


= ay & Pea L200 tO e= 360 


memea=l, (}=1); 


> Cod. + C3 See) tee T= 29 


eee taal = S(9) T° 242 51 


Moret -2, 3-0; 


c, = b, x d, + co) = cote t= 9360 

dg = a, xc, + dy =a 2eso0) + 0 = 720 
mor J=1; 

Cc, = by xd, + Cy —ecotget 120 = 17 i 

ds = a, x Co it dq, Seale eto Gi) == —) 1) 2 
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mee y=2, (j=i); 


C 3 = Da Oe d., tr C3 eel eet 9 = 1 0 

d, = a) X C3 iH d, =e) eto t= 7 
and d. remains unchanged, equal to l. 
Therefore; 

= LO 360 ayy alt 10] 

ieee) | 720 202 fou Heels (2-56) 
and the transfer function should be; 

2 
- SWS etl St LO 
2) ey a 
Peesst 702s +7 1stl . (2-57) 


Since the h.'s from the Cauer Second Form were used in the 
inversion algorithm from the Cauer Frist Form, the vectors 


C and D require reversing non-zero elements, resulting in: 


LO LO ey 360] 


Qo) 
il 


fee val yo LAC Ie (2-51) 


tO 
it 


and the correct transfer function is 


10s°+171s+360 (j= 


5°+715~°+702st720 


mcs) = 
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A digital computer program (FORTRAN IV) has been 
written for both Cauer First and Cauer Second Forms and is 


included as Appendix 3 with documentation. 
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ite LINEAR SYSTEM ORDER REDUCTION VIA PARTITION 
Ve UE SeAUERK SEe@OND FORM 


A control system,in general, can consist of many tens 
or hundreds of elements. In such cases, the problems facing 
the engineer include: (1) too many variables to efficiently 
handle; (2) the dimension of the system is too high to com- 
prehend; and (3) the modifications needed to meet required 
design characteristics are difficult to ascertain. A logi- 
cal approach is to seek procedural methods which reduce the 
order of the system to a manageable size yet maintain the 
basic characteristics of the full dimension model. 

A number of different methods for system simplification 
have been proposed for the reduction of high order dynamic 
systems to low order models of a more computationally or 
analytically tractable nature. The approaches used are 
quite different, but appear to fall into three main groups. 
The first is to ignore those modes of the original system 
which contribute little to the overall response. Davison 
[6] chose to neglect eigenvalues of the original system 
which are farthest from the origin, retaining only the do- 
minant eigenvalues and hence Fora time constants in the 
reduced model. The shortcoming of this technique is that 
many systems do not have any "dominant" roots [7]. 


Chidambara [8] essentially finds a reduced forcing function 
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so that the steady state values of the lower order model 
agree with those of the original system. The consequence of 
this method is that the approximate model give correct 
steady state values but incorrect time responses because 

the reduced forcing function does not excite the modes of 
the two systems in the same proportions [6]. Marshall [9] 
proposed the reduction of the state matrix by partitioning 
it and setting certain rate variables equal to zero in order 
to maintain the original steady state values. This techni- 
que, like Davison's, is based on dominant roots and, there- 
fore, exhibits the same shortcomings. 

The second approach is to search in some manner for the 
coefficients of a set of differential equations of specified 
order, the response of which is sufficiently close to that 
of the original system when both are driven by the same 
inputs. Sinha and Pille [10] proposed a reduction technique 
based on the iterative application of the matrix pseudo 
inverse algorithm [11] to determine a model of specified 
order which minimizes the sum of the squares of the errors 
between the responses of the original system and the reduced 
order model to a given input. The main drawback of this 
method is that the objective function to be minimized is 
restricted to be the sum of the squares of the errors. Sinha 
and Bereznai [12] presented a method which minimizes a 
specified error criterion for a given reduced order model 


of the original system, based on the pattern-search algorithm 
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of Hooke and Jeeves [13]. Although this method provides 
more flexibility than that of Sinha and Pille, it generally 
requires considerably more computational time due to the 
poor convergence properties of the pattern-search algorithm. 

The third category involves application of the theory 
of continued fractions. Methods involving this approach 
have been developed by Chen and Shieh [14]. 

Sinha and DeBruin [15] and Fellows et al [16] have 
established the fact that among the methods previously 
mentioned, the approach by continued fraction expansion is 


generally the best for linear model simplification. 


A. SIMPLIFYING A TRANSFER FUNCTION 

The general nature of a control system is that of a low 
pass filter. Therefore, model simplification should con- 
centrate on the steady state aspects of the response with 
the transient portion given secondary consideration. As 
previously shown in Chapter II, the Cauer Second Form 
exactly characterizes these miens. 


Given the nth order original system transfer function: 


ae i 
2 b;s 
a= 0 y 
sg -F 2 a sJ 
}=0 


where an mth order simplified model of the system (where 


mis strictly less than n) is desired, the polynomials in 


a7 





equation (3-1) are rewritten in ascending order; 


n-2 n-1 
; Bytb st...- -+++D, 48 tb 4S 
a 7 n-1,_n (ED) 
a ta, st..-- reed, 5S elas ts 
and expanded into a continued fraction: 
T(s) = ah 
hy rf S 
+ 
he S 
h. a S 
ae 
DOr 1 : 
hon (3-3) 


An mth order simplified model is obtained by keeping the 


first 2m quotients of the expansion, omitting the remainder; 


mis) = 1 


+ 
2m-1 


hoo (3-4) 


and performing the inversion of the truncated fraction. The 
inversion technique presented in Chapter II can be used for 


mits purpose. 
purp 38 





Consider the seventh-order system, representing the 
control system of the pitch rate of a supersonic transport 


pererart £ 10 J), described by its transfer function: 


So SUC OGSe 105 333) 


eee O00 SE 08333) 
5 /4+83.6355°+h0975 °+703428 148537035" 


+ 2814271s* + 3310875s+281250 eS) 


By continued fraction expansion: 


mos) = IL 
9.00036 + S 
=too2co t S 


-.036856 + S 


Wecticigey 2 Ula 2s S 


-00071478 
(3-6) 


Suppose a second-order simplified model is desired. Equa- 
tion (3-6) has fourteen quotients. For the desired system, 
the first four quotients are kept with all others discarded. 
The truncated continued fraction becomes: 
S2 00036 + S 
-.486286 + S (3-7) 


-.036856 + S 
5S dligd eS 5 


oo 





and converted into transfer function form; 


ecm eS 


Z : (3-8) 


T,¢s) = 
Ss t1.14644st+. 09941 


The block diagrams of equations (3-5) and (3-8) in the 

Cauer Second Form are shown in Figures 3.1 and 3.2, respec- 
tively. The unit step and impulse responses of the original 
and simplified systems are compared and shown in Figures 


eee and 3.4. 


Seo tale EQUATION SIMPLIFICATION 

The method of system simplification just presented is 
especially advantageous when converted into state space 
form. In Figure 2.3,a name for each state variable is 
given after each integrator, shown in Figure 3.5, from 
which the state equations and output equation can be 


directly written. 


ine) hy an See za 
Ppbeeny (hath, yh,(hjth,)  .--ho (h +h) Zip 
mises hy Ch, th,) ho(h,+h,+h,).--h, (h +h,th,) a 
Peace) Ba(hotheth,)...ho (hj th,t..ho 4) 142, 
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i 
il 
+ : ial and (3-9) 
i 
al 
M5 
“3 
Meeein, hy, he----ho_] : Cela 
a 
tL no 
Z meizet Br , and C = Lz . (= 1a) 


Simplification of the equations in (3-11) can be achieved 
Meebeartitioning of H, D and L, as indicated in Figure 3.6. 


The resulting mth order system becomes: 


Zz toes Dar, (=) 
where: 
eee ore =: nies 
.. een ch, tha) =... De Gaee tee | al 
Seen hth.) ».... Rom hy t+ +3 om | 
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ue 
1 
os = : 5 (3-14) 
a 
and Sp = as aor (3-15) 
254 
29 
aC = th, ny Sieh s howd : Gs—16) 
Z 
pm 


As an example, consider the seventh order system described 


Meeeene transfer function: 


C(s) _ 1441.535°+783195°+525286.1255+607693.25 
ey s/+112.045°+3755.925°+39736.735" See 


+ see S06 56s 775989". 195-+663656.255+617497.375 


Arranging the polynomials in ascending order and expanding 


into the continued fraction yields: 
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ets) 1 
R(s) 
ie Onsousi os + S 
ienOo li? ot S 


=O vom th S 


(3-18) 


2956560 .646 + S 


nNOOOOTS7 


and equations (3-9) and (3-10) are formulated from the 
quotients. A simplified model of second order is desired, 
therefore, the state and output equations are partitioned 
as indicated in Figure 3.7. The simplified transfer 


mUnetIOn iS: 


Mees) - -250367st1.035264 
aOCey —————-. .. ... i. (3-19) 
Seo U ST Oost .051966 


Unit step and impulse responses of the original seventh and 
Simplified second order systems are shown in Figures 3.8 and 
3.9 respectively. It is observed that the results of ex- 
pressing the seventh-order system by a second order model 
are satisfactory. 

It should be pointed out that a stable transfer func- 
tion may produce an unstable simplified function because the 
method of continued fraction expansion approximation does 


not necessarily guarantee a stable system. 
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Poles tGNn OF OPPIMAL LINEAR CONTROL SYSTEMS 
itn GpeeoChBBED EIGENVALUES 


A. INTRODUCTION 
Consider the control of a plant with dynamics described 
by a set of first order, time-invariant linear differential 


equations of the form 
x = Ax + Bu, (4-1) 


M@iere x is the nth-order state vector, A is the (nxn) 
plant matrix, u is the scaler control and B is the (nxl1) 


input matrix. The output is defined as 


Cx . (4-2) 


me Me 


<= 
iF 


where C is the (lxn) output matrix. 
A linear feedback control law is assumed, and of the form 
me = Gtx.” (4-3) 


i 


There are mainly two separate approaches in the deter- 
mination of the feedback control matrix, G*, corresponding 
to the system under consideration; 1 - optimal control and 
2 - modal. 

In the optimal control approach a performance index is 


considered which is to be minimized in the design of a 


fuulestartes are available or an observer or Kalman filter 
iS used to obtain the unknown states. 
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system. Assuming a performance index can be defined that 
represents most of the design requirements, "the solution 
to the optimal control problem can usually be obtained only 
by numerical methods that yield solutions to only a parti- 
cular problem" [17]. If solutions are sought to more than 
one numerical problem, simple performance indices must be 
defined, which often do not satisfy many of the design 
requirements. Therefore, the choice of a performance index 
must fall somewhere between a realistic criterion and one 
that is mathematically tractable. 

A quadratic performance index will be considered as a 


eriterion for designing linear systems, of the form 


J = 4 5 Cx" Qx + Ru-Jdt ; (4-4) 
O 
where Q is a diagonal non-negative definite (nxn) matrix, 
and R is a positive scalar. 
In the modal approach, Cs is chosen so that the closed 


loop system achieves the prescribed eigenvalues. Equations 


(1) and (3) together yield 

x = (A + BG*)x 

mM © and R are given in the optimal control approach, 
then the eigenvalues of the closed loop system are uniquely 


determined, which may not realize the required performance 


characteristics or desired degree of stability for the 
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system. Using the modal control approach, a feedback control 
matrix can be found that will give the system the desired 
eigenvalues. This matrix is usually not unique, and it is 
not possible in a practical sense to find one that is 
"better" than it's predecessors, since a performance measure 
is generally not known that corresponds to a given feedback 
control matrix. Therefore, it is necessary to find a method 
for determining the matrix, G*, that simultaneously satisfies 
the desired eigenvalues and minimizes a given performance 
index. 

In addressing this problem, Chang [18] and Tyler and 
Tuteric [19] have applied the root locus method to single- 
input, single-output and multivariable systems, respectively. 
The method lacks a rational computational procedure for 
determining the elements of the weighting matrix, Q, to meet 
a set of prescribed closed loop eigenvalues. Anderson and 
Moore [20] presented a restrictive method whereby a set of 
eigenvalues may be located to the left of a fine parallel 
to the imaginary axis in the complex plane. Chen and Shieh 
[14] presented a method using sensitivity analysis. 

Solheim's [17] method of a diagonalized (decoupled) system 
becomes complicated when the system contains either complex 
or multiple eigenvalues. Systems that cannot be diagonalized 
add further to the complication of the method. 

The method developed here takes advantage of the 
properties of the Cauer Second Form, is approached in a 


simplistic manner, and is easy to implement computationally. 
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Bee LRANSFORMATION TO PHASE VARIABLE FORM 


Consider an nth order linear system of the form 


e=SetfTf , 
eee Output 


d 


We 


(4-5) 


(4-6) 


Mmeverman [21], et al., have shown that if the system is 


controllable, then there exists a non-singular transformation 


matrix which takes an arbitrary state variable system to 


phase variable (canonical) form. 


where 





26 


e* ¢ @ @ 


(see Appendix A) 


(4-1) 


(4-2) 


ie eens = 0 
eee een 0 
Wee erene Bees ice. 
0 

SOG i 
a ear O 





If the system is not controllable, the phase variable 
(canonical) form may still be obtained from the system 


fmeeaaster function 





(4-7) 


Once the system 1s in phase variable form, Chen and Shieh 
[22] have shown that the equivalent system in Cauer Second 


Form is easily written as 
z= Hz + Vu (4-8) 


y= C*z , (4-9) 


~UllCl 


where the two forms are related by a linear nonsingular trans- 


formation matrix, P, 


ge= PR « (4-10) 


~ 


The matrix P is an (nxn) upper triangular matrix. 


oy 





The performance measure under consideration becomes 


CO 


oa = if ioe + mare alc C=) 
O 
where 
me <P>) Q PO 


From optimal control theory, the Hamiltonian 


lz Oz + u Rul + B'CHz + Vul , (4-12) 


mae Oe 


H 


where P is the set of Lagrange Multipliers (also called the 


o~s 


costate or adjoint vector). For the Hamiltonian to be 


globally minimized, assuming no bounds on admissible states 


and control values, it 1s necessary that dH/du = 0 and 
9*H/au->0. 

9H/du = Ru + VB = 0 (ae) 
implies 

u* = -R71y"5 mand (4-14) 

9°*H/au’ = R>0 , Gu=as) 
since R was defined as a positive scalar. Included in the 


necessary conditions for optimality are 


z= Hz + Vut (4-16) 


one Oe 


9H/az = -p = Qz + H C=) 


vt > 
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Combining equations (4-14), (4-16) and (4-17) yields a set of 
2n linear differential equations forming the canonical sys- 


mem in Cauer Second Form. 


Z H Naud, Z, 
s ~ T ; (4-18) 
P ace oe P 


It remains to be shown that this form is useful in obtaining 
optimal closed loop solutions that correspond to a set of 


prescribed eigenvalues. 


C. SIMILAR EIGENVALUES 


Consider the 2nth order cononical system in phase 


waptable form: 


Xx A -BR ~B x 
zs . (4-19) 
P ae “4 P 


It is known that this system possesses n eigenvalues with 
negative real parts and n eignevalues with positive real 
parte, and that they are located symmetrically about the 
imaginary axis in the complex plane [17]. The eigenvalues 


of the optimal feedback system 


x = (A + BGt)x (4-20) 


1 





are identical to those with negative real parts in the 
canonical system. It is possible, therefore, to focus 
on the canonical system where the dependence eigenvalues 
on the weighting matrices Q and R can be studied without 
solving the matrix Riccati equation. The problem is in 
determining a weighting matrix, Q, such that the system 
attains the prescribed set of eigenvalues. 

The canonical system in Cauer Second Form can be ob- 
tained from the phase variable form using the linear trans- 


formation 


Zeer Xs (4-10) 


where P is a nonsingular matrix. We have 
x A Sap e- x 
: . (4-19) 
P a a P 
: =. 7 
ir P 0 A aoe 
; = x (4-20) 
D 0 Ca leo eA 
ae 0 zs 
(4-21) 
0 oe 
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Let 


and 


ttt >« 


moe 


parvlatp! 
QP cp7tytatpt 
Eyre ; 
-H 8 
-BR tp! 
3 
_at 
-yoly! 
3 
_yl 
0 
3 
cp7lyt 
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vo > 


(4-22) 


(4-23) 


(4-24) 





where each sub-matrix of M, N, and P are known to be (nxn) 


~e 


square matrices. It is easily seen that P is non-singular, 
and that 
P 0 oe 
> a 
0 Goi) 0 pi 
ppt 0 
0 aie 
ii, 0 
= = I ° (4-25) 
0 af 


Therefore, 


me 


tO 2 
ero 
i 


M N (4-26) 


shows the similarity of the M and N matrices. Two similar 
matrices have the following properties: 


1. Their determinants are the same. 


det M = det N (4-27) 
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meee oneir traces are the same. 


eee tein g (4-28) 
k k 
a Mes = — Ns (4-29) 
3. Their characteristic equations are the same. 
detCLAI-M] = detfAI-N] = 0 , (4-30) 
where } is an arbitrary variable. Since their characteristic 


equations are the same, the eigenvalues of M and N must be 
identical. It 1s now known that the Cauer Second Form and 
phase variable system matrices are similar in that they 


possess identical eigenvalues. 


feeb VELOPMENT OF THE PRESCRIBED EIGENVALUE PROBREM 
Initially given is a known linear system with dynamics 

described by either a set of first order differential equa- 

frons Or its transfer function. It is desired to find the 


optimal feedback control, u*, such that the performance 


moasure 
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Pfilul eel (4-31) 


is minimized, where the eigenvalues of the optimal system 


are specified as 


r d r 


M2? 3? 


mrevaluation of the State (H) and Linear 
Transformation (P) Matrices 


Assume the transfer function of the known system is 


a (4-32) 





From equation (4-32), the Routh array is formed: 
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(4-33) 


c= [8 Bp Pe (4-33) 








d, d, Mee demer cite esi 3s'6: ari 0 
e, eo i 
0 
(4-34) 
Wied 
0 


In matrix notation, the Routh array becomes [rs .d, where 
Mea, 2, .2-, 2¢nti)) 


eae 2, ..-., ntl], 


and elements 


Bena k)+3, k  - 
mGniek)-4, kk 
Ment 2. J... nti] 


In general, the (2L+1)th row of the Routh array is the Lth 


row of Pp. 
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or 372 Se Sas Tl 
0 Dey Pen cette cece er 
0 0 Day crete eeees eo 
P = s 
“ (H=35) 
0 0 Oe Anca ce Rees 
Aries 


P being an (nxn) upper triangular matrix, will always have 
an inverse, eo. which can be quickly and efficiently deter- 


mined by digital computer. The H matrix formulation becomes: 


i> A p> (ices) 
Pci 2G: «<< ie 
eRe yy coe ci: ae 
a eee oe 
. Grae) 
d n-l,n 
een Danner a. ee Le. 


The elements of the H matrix can also be obtained more 


easily and directly from the first column of the Routh array. 
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Let 


ee 
ho = Giese) 
ee eee 


The h;'s correspond to the quotients of the continued frac- 


tion expansion of the transfer function in (32). 


ms) = 1 


(4-39) 


The H matrix then becomes as shown in Figure 4.1. The 
Geewlar pattern of the elements enable the H matrix to be 
obtained by inspection once the h,'s have been determined 
from either (4-34) and (4-38), or (4-39). 


The matrices V, and C*, are easily obtained: 


oH 


V=PBe : (4-40) 


f4 
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or 


wWIOJ PUODSaS Taneg UT 'H *XTAWeW 9IeIS “['b asANsTy 


(Eo cust. et Tyy2Cy 8 a 


(SytEytlyy ey e e e e 


(fytlyy Cy a «* e ee e 


Late 


U U Q e ® 8 e e e 


er ee ur reenact 


e (SyutEytlyy fy (futlyy uy Tycy’ 
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di 


mepcem = [h, h, h, ..-- hy] (4-41) 
B. Evaluation of 0 
g = (PT) Q Pm 

Sem 23 Co ate 

Bm 22 93 Con 

Q Q Osi. canoes Q 

ne. Ueee....... Can 


The canonical system in Cauer Second Form (4-23) has now 
been obtained, with the numerical values of the elements 
of Q still to be found. This system will be compared to 
a non-optimized system with the prescribed eigenvalues also 


in Cauer Second Form. The desired system in phase variable 


form is: 
x = Atx + Btu (4-43) 
yee = XxX (4-44) 
Formulation into Cauer Second Form yields 
x = H*z + VFu (4-45) 
y = E*z : (4-46) 


oo 





In a "nonoptimized" system (Q and R set equal to 0), the 


canonical system appears as: 


Z, H* 0 Z 
= ; Gea 7 ) 

pala 

Dp O - (H* ) Pp 

where 

H* 0 
= NZ (4-48) 

0 -cH#)! 


possesses the n prescribed eigenvalues with negative real 
parts and n eigenvalues with positive real parts, symmetric 
about the imaginary axils in the complex plane. 

By equating the characteristic polynomials of N and N¥, 
(4-23) and (4-47) respectively, it is now possible to 


determine the elements of the Q weighting matrix. 


det[sI-N] = det[sI-N*] (4-49) 
er vRotyt 
det ; = 
Q Pye eie 
Sioa 0 
det 
0 SG) (ies0) 
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PDE LERMINATION OF THE OPTIMAL CONTROL LAW 

Finding the elements of the weighting matrix, Q, ales; 
obviously a non-trivial matter for all but the lowest order 
systems. The method developed for obtaining these values 
is based upon a succession of matrix building blocks, which 


are individually computationally simple. 


Starting with the matrix, H, 


lS ths] (4-51) 


define a new matrix, T, 


where 
f= f.. -— h.. eee (4-53) 
i 53 : 
and 
n 
> ae : rterees ; ieee it (4-54) 
Tal j=itl a7] 


The matrix, T, therefore, is an (nxn) upper right triangular 
matrix, where the diagonal elements are equal to the sum 
of all other elements in the same row. The next "building 


block" matrix, G, is defined by: 
G = {g..]. (4-55) 


= 1.0 (4-56) 


ae 





Memmi il, 2, ...-, nj} , 


aa OC }-2, 3-1 (4-57) 
Setgeeyelt, 2, ..-., n-l] , 
n-j+2 
ae. = z CE ) (4-58) 


. x ° 
Tag) eae os ek ,j-1 


mummemeic, +,..., x] and iE[1, 2,..., n-j+tl], where the 
index } is held fixed for each summation over the index i. 
The matrix, G, is an (nxn) upper left triangular matrix 
characterized by the first column being all ones. One more 
matrix needs to be defined at this point. Let the matrix, 


W, be defined as: 


Bee Paid (4-59) 
where 

; oe ; 

Mek 7 St X84 (4-60) 


fer i, |} and k e[1, ey aie oe 


The matrix, ie is therefore a tridimensional array with each 
"level" an upper left triangular matrix. Examples of the 
matrices, T; G and W are shown below. Although not evident 
at this point, the G and W matrices will be used heavily in 


obtaining the values of the elements of the Q matrix. 


V2 





From the linear transformation 


foes >) QP =, Geen 


it is observed that once the element Qs5 is known, the 
remaining elements in the same row (column) can then be 


obtained through a process of 


Fe ze c 


11 12 1 i 
0 vo t 53 Reeua eran eines ton 
0 t 33 eee een tah 
0 
eo = 
(4-62) 
ce 0 CR octiioeemee ater et gus Tan | 
il ThA a oo 2a 
il to59 ee oOo 0 
1 t 33 er 0) 0 
= 
En-2,3 (4-63) 
i en: 0 0 
Hl 0 0 0 0 


nS 





(h9-h) 


eee? @ 


us 
ag 
u 

=£ 

C qt 
fa 3 
U 

-C 
eel! 
Ca 
u 

Ct Gat 


oy 
cw 


= ee 
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linear combinations of previously determined values. The 


Q matrix is thus found in the following manner: 





a Q11 =. 17 .3 213 O11 eT eer: Qin 
Oop i Qo 9 7 + Qo se on 
O33 on O34 es an 
an = 5 lun 
ng 
(4-65) 


where for an nth order system (n>1) 


2h 2a 
0 _ ae il 
a 2n 
Te Ge.) 
es 
i al 
= nm (;*)° [(HH,)?=(hoh,)°] * (4-66) 
aig ale 


+ E - z 2 2 
fOr n=1, Q44 = (HAH, ) = Chahy) : 


Vs 





~ = a + 
ml Ee ik (4-67) 
4 ee pee 
mem if] or jfn , 
eG n-l ee 
©. = > Q.. (4-67) 
in Peal 5 
and 
Ee n ote n n 
ite, 6) =] J =a jen JJ 
n n n n 
meek | Oy Boge (H..-H..~) - £ Bea Cee ) CG ) 
izl Fel -~ll JJ LJ L=l j=] LJ a3 1) 


The values of the diagonal elements of 0 Cother than Q54) 
generally involve varying linear combinations of already 
determined values, more easily expressed in terms of the G 
and W matrices, rather than the P and H matrices. To aid 
in determining the values along the diagonal the following 


labels are provided for G and W. 
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S S S 
Ik S11 815 813 ae ee 
2 Sor S99 S53 EM ene 
Memeet)=— 2320S 833 te 
4 
ee) on-1,2 
Enl 0 
-N-1 N-2 
Ps 
i Pe tan oe 
ee ihe ie 
Q32 213 Wood 
43 ane Weoq 
n-1 ei. 4 et. 
en Peeler t 





77 
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(4-70) 


| 


(4-71) 





These two arrays will provide the coefficients of each ae 
sh" necessary. Assume, for example, that a 5th order 
system is being considered in equation (4-50). The results 


of the expansion of both sides of the equation results in: 


n : n . 
ao fa ey = bY 7 os b.s** (4-72) 
1=1 L=1 
where a, = ells bo = 1. By equating as 4 and Ds 4 core 
meme. ..-., ni, it is possible to obtain Qaas 
a and b. are the coefficients of ee From (4-70), it is 


observed that the only element in the 3 column that is 
non-zero is Eis: which appears in row 1. The 1 "indicates" 
that it is necessary to only look in "level" 1 of W, in the 


column corresponding to ae This yields only the single 


0 


element Wisi: Therefore, the coefficient of Q318 will be 


Bee < Wie, (HEF > 


Wheat remains to be determined are the coefficients of =o 


mot involving the veggie Because of the symmetry of the 
eigenvalues of both the system being "optimized" and the 
system with the prescribed eigenvalues, the remaining 


coefficients (those not involving a aa are easily 


determined from: 


7s 





detCsI-H] x det[sI H’] 


Ce 755 
and 
det{sI-H*] x det[sI+(H*)"], (4-76) 
which give: 
n : n 
Meee ss )( £ (-1)%a, s*) (4-77) 
i=0 fei = 
and 
n : n 
( £ a,#s7)( £ (-1)%a,#s*) (4-78) 
= P= 
mimere K=i+] for n even and k=1 for n odd. The indicated 
multiplication in (4-77) and (4-78) result in: 
n : 
meec-1) oes-2 (4-79) 
i 1 
1=0 
and 
n ; 
gS  (-1)*g,%5%? (4-80) 
: a 
1=0 


respectively, where the same conditions are imposed on "k". 


Returning to the 5th order example, the coefficients of 


Ss are -o, and -o,*. To obtain Q44 is now a matter of 


solving the equation: 


~ 


(B35 ° Wisy)? Q3 7 [9 = 79Q" 


0 (4-81) 


1S 





ote 


Ter 
Q ra (le ae (4-87) 
Jl (845 * 835) 


With Q14 known, it 1S a simple matter to obtain the 
remaining elements in the first row (and column) using 
(4-67) and (4-68). 
ikon rind a find all the coefficients of ae 5° results 


in three separate ways: 


@> s xs (4-83) 


The multiplicand indicates which columns of G, (4-70), is 
of interest. Any non-zero element in G tells which level 
of Wis of interest. The multiplier 1s the indicator for 


the column of interest in the array, W. Therefore, for 


B15 * “Wy 37 4 tye ae ee sacm Q,3) (4-84) 
for st x st. 

eee 74 21 otal Qi») 

Bou * Win Qo * Wono Qo 9) (4-85) 


80 





for Ee xX Ss 


fe 


x Wis, 47? 


$153 


x CW oe 


82,3 152 201? 


pee “isa 231? - caer 


~~. 


Obtaining Q50 1S now a matter of solving the equation given 


by the combination of (4-79), (4-80), (-84), (4-85) and 


(4-86): 
3 ~ Z 2 s 
: : . + ; Keer 
E15 . Mee ae a Soe ee Sea 
3 Bs 
for Q59° 


Once Qo has been found, the remainder of the 2nd row 
(and column) can be obtained using (4-67) and (4-68). 
In a like manner, the respective equations to be solved 


nor Q33> Quy and Qe are: 


ale 5 . 2 ; ; 
Be: ee 7 Wa 3 
: f 4 2 a 
‘ ‘gs °453 i213 854 ‘ ey? pees S54 
° if : 
; ae P4551 asi Ea a (4-88) 





8 5 ~ + uy Rs 
—)6|CUS Eee Ones + € 2 og Dee ree Oe a ) 
$=1 Woe ea aa ao S= ae seal de ae 
S 3 ~ 
mcen 86 Tee Or Oo. = oF , (4-89) 
j=l 1a ea egy yee 3 3 
° 5 = 
ie ¢. Meer Oe )r— 0, =-—0, = ... (4-90) 
321 Bispl eal lets) a 4 nm 


The underscore in (4-88), (4-89) and (4-90) indicates the 
term that contains Q39> Quy and Qos respectively. 
The entire Q matrix 1s now known. By the inverse 


transformation of (4-42), 


a= P 


~s 


(4-91) 


2D 2 
t'U 


eeeeong with A, B and C of (4-1) and (4-2) are used to 


obtain the matrix Riccati equation steady state solution. 


0 = KA + AK - KBR pix +0 . (4-92) 


~~ ~e 


Once K has been determined, the optimal control law is 


~~ 


given by 


i = -R7 BKx a : (4-93) 


~~ fe we 


Once u* has been determined, an inverse non-singular trans- 
formation can be performed to take the phase variable form 


r fais now known that Ea aaK in (4-20). 


or 





back into state variable form. (See Appendix A) 


mee LLUSTRATE EXAMPLES 


imeecda Order System (Third) 


Given the system 


xy 0 i 0 x 0 
Xo = 0 0 ay Xo + 0 u (4-94) 
we 3 -13 -19 5 th X4 ap 


(4-95) 


find the optimal control law, u*, such that the quadratic 
performance index (4-31) is minimized, where the eigenvalues 


of the system are specified as 
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Forming the Routh array yields 


13 ine 7 i: 
1 0 0 

iS i Ih 

-7/19 -1/19 0 

20/7 a 

ty 3-0 0 

1 

0 


from which the third, fifth and seventh row are extracted to 


form P. 

19 iy 1 
= 0 S07 a: : 

0 0 a 

and 
1/19 ~49/570 1/30 

Ea = 0 7/30 -7/30 

0 0 al: 


From (4-97) and (4-98), H and V are calculated 


~13/19 637/570 ~13/30 
H=pAp7t= | -13/19 -63/19 9/7 
-13/19 -63/19 26 
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(4-97) 


(4-98) 


(4-99) 





Y = PB = i (4-100) 


xy 0 all 0 x 0 
Rp = | 0 0 il Ps, + 0 u Gs= One 
we -72 -54 = VS Jf 


and 
y = Ll 0 0] xy (4-102) 
2 
cory 
From the Routh array 
eZ 54 ie a 
a 0 0 0 
54 iS 1 
mals / 54 -1/54 0 
Hee / 13 1 
115 0 
i (4-103) 
0 


So 





extracting rows 3, 5 and 7 yields the P* matrix: 


54 ie 
pe = 0 ive ns 
0 0 


from which H*' and V' are calculated, 


~4/3 
H¥=P#as(p#)~* =] -4/3 
any 

Se 

i 

i 


Bios S 
-2574/621 


-2574/621 


1/54 
0 37 ilo 


0 0 


Sey aS 
1980/1495 


-11254/1495 


Formulation of equation (4-50) yields 


86 


SiH 


-169/6210 


Wes 
Sos Lis 
IL 


CEST GH) 


(4-105) 


(4-106) 





pumeoy 19 #«8-637/570 13/30 1 il 1 


3/719 s+63/19 -9/7 i li iL 
iy 19 63/19 st+3 i lt il 
det Qn4 Q35 Q33 s-13/19 -13/19 -13/19 
Qi Qo Q43 637/570 s-63/19 -63/19 
3 Or one eller mS) 7 S28 
Bee/3s —-676/345 FO MALS 0 0 0 
4/3 Ss+2574/621 -1980/1495 0 0 0 
3 Dee 6 2) s+11245/1495 0 0 0 
mock |G 0 0 Be ea ay) a et ys 
0 0 0 676/345 o&2574/ -2574/ 
aa S2ZE 
8 0 0 -72/115 1980/ s-11245/ 
1495 1495 | 
(4-107) 
It 18 desired to determine the values of Q. Brute force 


~s 


enumeration of the determinants and equating coefficients 
of like powers of s would eventually lead to the solution. 
Using instead the method developed, from H(not H*) evolve 


the T, G, and W matrices. 


PaoOGG #24333 2.5667 
T = 0 He2eas? 4.42857 


0 0 0 (4-108) 


87 





< 


i ye Ou00 


G = az eco Of 
I 0 
ae 7.0000 
and W = -1 ee OO] 
-~1 0 


Now add the appropriate 


for G: 
Z a 
S S 
1 ak 72 
2 a}: Lae OL 
3 i 0 
for W 
Z 1 
S S 
Qe4 -J] Po 
Qs5 -1 “u.2857 
Qi3 L-1 0 


Still necessary are the 


n 
gp 


° n 
a 
a-s )€ F 
zl! - Le 


0 (0) 


k a 
(-1) cher 


19.0000 
Q 
(4-109) 
0 
-19.0000 
0 
(4-110) 
0 
labels as in (4-70) and (4-71) 
0 
S 
19.0 
0 
0 C4Y-J11) 
0 
. 
-19.0 
0 
(4-112) 
0 


coefficients in (4-79) and (4-80): 


y= (-1)*o,s7* (4-113) 


im a 


1=0 
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Bemen odd (n=3), k=itl, 


n : nN ° 5 Ti - 
+ 
(eee a.s”)( a (is lyst) = ye} (-1)*"*a, 21 (4-114) 
1=0 Sa = 

Ws -4195+13)(s°-75-+19s-13) 

= =o = iene te 179s° - 169, (4-115) 

and 
3 . 3 : : 3 ‘ 
%, + f, f, 
( E a,#s*)( £ (-1)7*ta,#s*) = 2 (-1)7*t0, #577 (4-116) 
n=) 310 nei) 

fe les -+545+72)(s°-13s-t+54s-72) 

2 BP = ee + LO4uus? - 5184 . (4-117) 


Now, by equating coefficients of ay it is possible to 
obtain ..- From G, the only non-zero element in the ge 
column" is 19, which corresponds to row 1, and therefore 
level 1 of W. From SC uel iy W, the only non-zero element 
in the 169 column" is -19, with row label Q01° The 
coefficient of Q1, as obtained from (4-81) is 19(-19) = 

Soot. The solution for Qs 4 is obtained from equations (4-81) 


and (4-82): 
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-361 1 - 169 = =-5184 


Q34 ood £ 169) / -361] 


bee 


MeO. . s (4-118) 


~~ 


From (4-65), (4-67) and (4-68), Q15 and Q33 are successively 


obtained: 


“1 


x z . By ole 


- -22.69 (4-119) 
Zee 


Wi 


Seis oie = C= 227.69 ) Caw OS. (4-120) 


oo 1 ~ 12 


The next power of s which results from expansion of (4-107) 


1s 5*. Equating coefficients of s*, it is now possible to 


obtain Q59° a results from the products s xs? Sx. 


and s?xs?, Therefore, from the development starting at 


(4-83), the equation to be solved for Qo9 is: 


1 3 ne 2 2 = 
er ¢ PRN oon.) © (2 ig. esos.) 
j=1 3,3 izl ely jt =n 332 izl 12] jt 
g i z 
r pepo.) (Cl UW eee Cr (4-121) 
j,i ea 13] Q54 i iL 
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Substituting known values, equation (4-121) becomes: 


213¥11%) aaa Qe ene - 
852%) 91 a UO as * B50 W102 Qo * Wop9 Qo) 
ead "131 Qay? as 54) * Bays 85? 
ao, = o,* (4-122) 
Meme =1(13.892) = 1(-22.69) - 1(8.798)] 
feo £7.0(13.892) + 4.2857(-22.69)] 
mees7 (7.0(-22.69) + 4.2857 20 + 1£-19(13.892)] 
Pee —19(-22.69)] + 10-19¢8.798)] + 179 = 1044 
Qo 5 = 84.155 (4-123) 
From (4-68) 


one 


ee seo 222 ~ ~212 ~ Yoo 


meee 2 69) —- (84.155) = -61.465 (4-124) 


gal: 








e is the next power of s obtained in the expansion of (4-107). 
Equating coefficients of s*, it is now possible to obtain the 


equation from which Q33 can be found. s" results from multi- 


plying as. Therefore, from G and W: 


Prpe-1(Q,,) - 10,5) - 1Q,,)] 


+ 


Boz b-16Q,4) = 1(Q,5) = 1(Q,3) J 


+ 


G3,1-1(Q.,,) - 1(Q35) = 1(Q33) 


i -5,* 5 (4-125) 


where g,1> 85, 4nd g3) all equal one. Therefore, (4-125) 


becomes 


(ew O65 = -Cn (4-126) 


t 
1h &3 
(8 (4 G 


isl jj 


Beiving (4-126) for 033° 


‘. = -a,% + 0, - z eS hes (4-127) 
4 5 
since 
5 ee 
Qin = > FE, bt (4-68) 


OZ 





Q.. 20; (4-128) 


and (4-127) simplifies to 


: Pee 
fee 2 ~ Fp - % Qa, (4-129) 
k=1 
Pee i-8.798+61.465 = 102.667 (4-130) 
Q is now entirely known. 
13.892 Oe GG 8.798 
Q =| -22.69 Saag ZBI S Geis 
8.798 Bee tnES lO eee 
From (4-91): 
5015 0 0 
Q=P QP = | 0 865 0 (4-132) 
. 0 50 


To find the optimal feedback control law, u*, it is necessary 
Memeolve the matrix Riccati equation (4-92), where Q, A, B 
and C are as given in (4-131), (4-94) and (4-95). The 


solution yields: 


a3 





u* = -£59 35 6 J xy (ue 33) 
0 
= 


or u* = -59 a 315 Xo 6 X4 


With G* as given in (4-133), x = (A + BG*)x 


becomes 
xy 0 IL 0 Xa 
ms - 0 0 i 7 (4-134) 
2 2 
; -72 -54 -13 
ea og 


Comparing (4-134) with (4-101), it is observed that (AtBG*)é= 
A*; therefore, the desired system has been realized with the 
set of prescribed eigenvalues. 

2. Even Order Example (Fourth) 

Consider the linear system represented by the transfer 
monet 1lon: 
Zz 
| (3) a (4.135) 


5 t+125°+h85~°+805th8 


the system eigenvalues are -2,-2,-2, and -6. 
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This system could be expressed in phase variable form (4-31), 
thereby obtaining the transformation matrix, Ee the system 
can be realized in Cauer Second Form (4-36). Instead, in 
this example, it is desired to obtain the quotients which 
result from partial fraction expansion of the transfer func- 


tion. In Cauer Second Form, T(s) becomes: 


eee 
3. 4 


ines) = =) 
PmoaogSsGtiesutl 2s +s 


(4-136) 


2 





where the quotients, n> are: 


hy ewe. 41176 he Sok 

h, oO O245 he = .146914 
h weet 620 7 ho = ~781.808 
hy = ~,627918 Ne = =, 924 ONT 


substituting these values into the matrix in Figure 4.1 


yields 
e378 —.88647 20741 =. 0S422 
ms78 2.01997 -.47261 Ue oo 
— . cealeg. 
c mrs72o 2.01997 ZG ae ~.48125 


meses 2.01997 ZO) @ 6.35004 


Again, V is (and in all cases considered will be) a vector 


t 


of all ones, 


(4-138) 


i<j 
i 
PrP FP 


It is desired, as before, to find the optimal control law, 
u*, such that the quadratic performance index (4-31) is 
Minimized, and that the optimal system realizes a set of 


prescribed eigenvalues. The eigenvalues are given as: 
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aS ae ae (4-139) 


With the zeroes of the transfer function (4-136) the same, 
the transfer function of the desired system with the pre- 


scribed eigenvalues becomes: 


2 
T,#(s) = —-——S—2S" st (4-140) 
eergsetL38s t435st450 


By continued fraction expansion of T*(s) the quotients 


obtained are: 


H, = 13.2353 H, = -1077.43 
H, = .107635 H, = .004834 
H, = -71.3206 H, = 396.926 
H, = -.088175 H, = -.024294 , 
meom which H* is: 
1.42458  -1.16702 .06399 -.32154 
ieese 8 5,121168 =O ee ee sell 
H < _—- = 
1.42458 5.12168 -5.48975 27.58655 
1.4258 5.12168 -5.48975 17.94349 | (4-141) 


oF 





Formulation of equation (4-50) yields in matrix notation: 


a 
Spee L. sl-H* 0 


A, 
(D 
ct 
iN 
ul 
QO. 
Ay 
ct 


Q sI+H 0 See 
(4-142) 
which is to be solved for Q. 
As before, we formulate the T, G and W matrices 


according to equations (4-53) and (4-54); (4-56), (4-57) 


and (4-58); and (4-60), respectively. 
iiais yields: 


me? 00 2.9064 2.7093 6.3500 


0 Geeot’ 3.3893 -'6.2721 
T = C=) 
0 0 6.8313 6.8313 
| 0 0 0 0 
Pal 12.0000 46.5882 67.2941 
1 9.6614 23.1534 0 
= 
: 1 peo t: 0 0 CH-14u) 
1 0 0 0 


a5 e e e@ e e 
mme notation L indicates an (nxn) matrix with all elements 
Mnaity . ‘ 
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-1 ZO 00 0 -46.5882 67.2941 
-1 sc) leo) aig =23,,.1534 0 
W = (4-145) 
-1 Sr cones 0 0 = 
pea 0 0 0 4 
With the proper labels, (4-70) and (4-71), 
for G: 
5° ae a ee 
il a: 12.0000 46.5882 of 5290: 
2 i Wao o ls 23353 4 0 
(4-146) 
5 el Gos 13 0 0 
eis 1 0 0 0 | 
for W: 
3 2 1 0 
S S S S 
ma & 
S ¢ 
Q51 -l 12.0000 46.5882 67.2941 
oF -l 9.6613 23.1534 0 
a (4-147) 
Q; 2 -1 Books 0 0 
7 
Fi, Si ane 0 0 
G and W contain all of the coefficients of Q.58°" Resuie Gane 


from the expansion of (4-142). Still to be found are the 


og 





coefficients of s 


2k - 


Therefore, from (4-79) and (4-80): 


and 


We les F485 °4+80sth8)(s'a125°+h85--805t48) 


Me es” — 480s’ - 17928* + 2304, 
O54 = 2304 
om = -48 
o, = -1792 
ili S 
Oy = 1 
0. = 480 


fr 19s°+138s-+435st4S0) x 


3 


Wee 9s leis 4 35ee550) 


6 2 


= eos - 3414s" = 690205) © 202500» "<. 
o,* = 202500 

a -85 
oi -650295 ; 

Cia 
ill 3414 


100 


involvi ~~ na es. 
mei) lving any Qi? the g. 8s and a," 's 


(4-148) 


(4-149) 


(4-150) 


(4-151) 





Determination of the oi results from equating 
coefficients of like powers of s from the expansion of the 


determinants in equation (4-142). The Seave are obtained 


in a successive method (4-65) by starting with the aC 


coefficients, then continuing by equating coefficients of 


<a kell, an —1 i), in seandent fashion. 


Bor 5°; 
from G, the only non-zero element is B14? corresponding to 


row 1, therefore, level 1 in W. The only non-zero element in 


level 1 of W in the ey column is Wout: Equating coefficients 
yields: 

(Zi4 x Wiuy? OTe 8 o,* : (4-152) 
polving: 

~  _ _202500-2304 ~ - 

MMMCGTEDOGT) (67.2051) ~ '172080e ae 
and from (4-67) and (4-68), 

Qin = Ay = 788. 95327 

Q13 = 24, = 48.14819 

a & (4-154) 

Q, = %, 7 -3-40294 


neal: 





for coefficients of =o 


Cena ea : (4-155) 
Bolving: 


oe = 296.94152 (4-156) 


and from (4-67) and (4-68), 


=26or. 10162 


Q53 = 232 


Qoy Quo on Loa? (ie 757) 
morecoefficients of a 


2 Ly ~ 3 
Lb g. a W. 
fea? 4=2 


2 
meee. «5 W.5. O-.- + o, = o,*% (4-158) 
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polving: 


eee= 252.24159 
end from (4-68), 


Q., = 5 = -135.68816 


. 6 
mor coetficients of s : 


4 4 ~ 

Mewes W.-. Q.. + 6, = O,% 
321 es =) se 5 3 
‘ uo 

ee + cae 
Quy = Paap cl nes 
i+ j#8 

~ 3 
es eae . + Cae 

Quy wey 23 oe 

OL eel 0.3/7 73 


we 


Q is now entirely known. 


Gis 


(4-159) 


(4-160) 


(4-161) 


Cia 2) 


Cu 6 3) 


(4-164) 





Ph. 20803 mieioe 5327 48.14819 = egymeh 


9 = | ~88.95327 296.94152 -263,70162 55.71337 
; 48.14819 -263.70162 351.24159 ~135.68816 
-3.40294  55.71337 -135.68816 120.37773 | 
(4-165) 
From equation (4-91): 
200196 0 0 0 
r 0 63233 0 0 
= ae Geto) 
0 0 2934 O 
0 0 0 3h 


Substituting matrices Q, R, A, B and C into the Riccati 


equation, and solving, yields the optimal feedback control 


law, u*; 
u* = -[402 S10 J 7 x4 
oy, 
X. 5 (4-167) 
ah 
where G* = -[402 300 90 gale (4-168) 


The optimal closed loop system, : = (A+BG* )x, 


in phase variable form is: 


104 





| xy 0 1 0 0 x 
Xo 0 0 al 0 Xo 

= (4-169) 
X4 0 0 0 I. X. 
| Xy -—450 +435 +138 -19 xy 


which realizes the desired set of eigenvalues. 
3. Higher Order Example (Seventh) 

The previous examples represent an odd and an even 
ordered system, illustrating the minor differences in com- 
putational procedures. It is observed that for low order 
systems, the calculations can be done, relatively easily by 
hand. Higher order systems require, laborious and tedious 
computations. Appendix C provides a digital computer program 
which yields the weighting matrix, Q; with the only required 
input being the transfer functions of the known and desired 
eigenvalue systems. 

The following seventh order example utilizes the 
results from the program given in Appendix C. 


Consider the linear system given by its transfer 


mane tion: 
Ta(s) = ———,;@—__*2 5 (4+170) 
seroeus t40.4s +11678sS F233-0s 323.25 


je OO OG ety 20 310 
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It is known that the characteristics of the desired system 


ere such that its resulting transfer function is given by: 


Zt io! 6.6 


7 6 9 Ly 


To*(s) = 
s +15.4s +101.64s° +372.68s 


+819.896s° 


(4-171) 


+1082.26272s°+793.659328s+249 4357888 


The diagonal Q matrix was determined from the computer pro- 


~~ 


gram yielding: 


Q,, = 45834.21273428 

Qoy = 89780.21841734 

Q@,, = 55970.39786199 

y= 19170.7157376 

Ome = 32959.23664 

Qeg = 494.9776 ee 
ee = 33.68 


eek and the state and output matrices representing the 
transfer function in equation (4-170) in phase variable form 
were substituted into the matrix Riccati equation. The opti- 
mal feedback control law, u*, resulting from solution of the 


Riccati equation is: 
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Meee —127)1 .435789 xy =o US oo8 sg Xn 


mod .06269 X= Seo 2I096 Xm Zosnoo9 74 Xe 


oom. 2399914 Xe = 6.39999944 Xa 9 Cis) 


where the _ matrix G* is: 


ye 8 121.435789 $05.649319 759.06269 586.29596 


255.879974 61.2399914 6.39999944]. (4-174) 


The optimal closed loop system, 


oS s 
i 


(A + BG*) x 


cee cx 


expressed as a transfer function is: 


Y(s) _ Dinou eis ae 
3 ( ) ce 5 
ae 5 /+15.39999974s°+101.6399914s” + 
4 3 
SOG 97 Nien tec Scat 


1082.26269s°+793.659319s+249.435789, 


which realizes the desired transfer function with the pre- 
scribed eigenvalues. 

Figures 4.2 through 4.13 show the impulse and step res- 
ponses of the three previous examples, both before and 


after compensation. 
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ve CONCLUSIONS AND DISCUSSION 


A method has been shown for determination of the state 
weighting matrix in order to satisfy a prescribed set of 
eigenvalues through phase variable state feedback. From a 
strictly mathematical viewpoint, this technique requires only 
a knowledge of matrix algebra. Every attempt has been made 
to avoid the necessity of inverting a matrix. The intro- 
duction of Chapter IV made known the fact that previous devel- 
opments inthis area have suffered the main drawback of res- 
triction. The author believes the method presented here, 
using Cauer Second Form, overcomes many of these restrictions. 
It presents a rational computational procedure for determina- 
tion of the weighting matrix, Q; the system elgenvalues are 
only required to be in the left half of the complex plane as 
opposed to the left of a line parallel to the imaginary axis; 
and the method is no more complicated for multiple or complex 
eigenvalues than a system with linearly independent eigen- 
values (or eigenvectors). 

It should be noted that some authors "define" the 
eigenvalue(s) of a matrix to be only the real root(s) of the 
characteristic equation. In the development that has pre- 
ceded in this thesis, all roots are considered eigenvalues 


we tne associated matrix. 
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The algorithm derived in Chapter IV, is designed as a 
basis for future research. In particular, if the designer 
is working with nth order systems where n is relatively 
large, it may be advantageous to look at using mth order 
Simplified models (m<n) of each system by a partitioning 
scheme similar to that in Chapter III. Again, it is em- 
Phasized that reduced order models do not necessarily yield 
stable systems. If the simplified system retains the basic 
characteristics of the original system, especially in steady 
state, then this would appear to be a reasonable approach. 

A parallel approach could also be inventigated regarding 
multi-input multi-output (MIMO) systems, treating each 
element of the system transfer matrix as an individual 
ePanster function. 

To the author's knowledge, no work has been done in the 
digital or sampled data areas involving continued fraction 
theory. This area should be considered due to the increasing 
use and need for digital control systems. 

These topics represent just a few of the areas available 


for future work. 
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APPENDIX A 
Consider the system 
mee Ax + Bu, com) 


where x is an n-dimensional state vector, u is the input 
manetion, and A and B are time-invariant (nxn) and (nxl) 
matrices, respectively. The phase variable (canonical) 
system representation is defined as 


yo = Av + Bu, Ce 


where v He an n-dimensional state vector and 


0 ab es cet: 0 0 
0 0 1A eee 0 0 
0 0 0 0 0 
A = 5 B = (A=-3) 
0 0 (egos cetes ones AL 
Ay Wan “Agee reese wae, 2 


The systems represented in (A-1) and (A-2) are said to be 
equivalent if and only if there exists a non-singular matrix, 


K, such that 


~~ 


ie KV (A-4) 


~w 


te 22 


—_ == 


7 





Kalman [23] has shown that a necessary and sufficient condi- 
tion for such an equivalence to exist is that the system 
in (A-1) be completely controllable. 


The controllability matrix of system (A-1) is defined by 


memes | AB: ACB; .... 2 Anm4B) (A-5) 


~~ ~~ 


or in an equivalent manner 


ime [Le e 


1 > el, (A-6) 


meeme the (nxl) vector e. is recursively defined as 


@ Ae. 


ey 7 Ae. > 2, = i. (A-7) 


~ 1 


aA 


The controllability matrix of system (A-2), E, is defined 
in a similar manner with A and B. Since there is only one 
control input, a necessary and sufficient condition for 
Semepollability is that the (nxn) matrix E (or E) have an 
inverse. 

Silverman [20] has shown that if the system in (A-1) 
is controllable, then the transformation matrix, K, 1s 


determined by 
[enn rE, (A-8) 


where 


es 








a ay, eT cee an 
ay a. er eee Ree iu 
, a. ae Me 2. eS 2 0 
ae = 3 (A-9) 
ile ME Poe oes. feo sc 0 
0 ae ee 0 
ea 
2 
23 
and a= ‘ = -E* AY B. (A-10) 
An-1 
a 
n 


im@emelements of a are the coefficients of the characteristic 


polynomial: 


n i-l 


Getpsou—-Al = det{SI-A] =S + £ a.S 


t=1 + 


(A-11) 


The matrix inversion in equation (A-10) can be avoided by 
using the Leverrien-Fadeev method for calculating the 
coefficients of the characteristic polynomial. Once the 


aw 


coefficients are known, a 1s written by inspection. 
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Rane [24] presented a simplified procedure for finding 
the transformation matrix, K, requiring no matrix inversions. 


Substituting equation (A-4) into (A-1) and premultiplying 


equation (A-2) results in: 


x = AKv + Bu Cle?) 


~s ~w ne me ~e 


=exAy + KBu  . (A-13) 


ne Ne ee 


Comparison of equations (A-12) and (A-13) yields 


wn 


AK = KA (A-14) 


me me 


and 


Be-> KB. (A-15) 


~~ me 


Partition K into n column vectors, each (nxl), so that 


Le = Tk, Ko wee ki: (A-16) 


Substitution of (A-3) and (A-16) into (A-14) and (A-15) gives 


Alky Kp Kg v+++ KA] > 


0 lees, 0) 0 Ce) 
0 0 0 0 
[x Ky k3----- KJ 
0 Not eae 0 uh 
~a, “Ao-> ae =e 
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and 


Bs [ky Ky K one kJ 0 = ies : (A-18) 


From (A-17) and (A-18) 


i 
> 
4 

ey 


al a Sa are i 


enim 2 i 


Ky = eh Ko + kK Bo > (A-19) 
or, in general, 

oem ee en-itl Sn Sn-itl a 
Megeietl, 2, ..-., n~l]. The column vectors Ki> ee ia are 


found in a simple recursive manner and completely determine 


the transformation matrix. 
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APPENDIX B 


INVERSION OF CAUER I AND CAUER II FORMS 


This program was written in FORTRAN IV, and requires 
minimal input. The only information required is: 

1. the order of the system 

2. which inversion is required 

3. the quotients from the continued fraction expansion. 


Multiple data sets are possible, and input in the following 


format: 
Card Columns Description Format 
it 1-3 M = the desired order transfer 13 
funee der 
2 1-20 h,, the first quotient of either D7 Oes 
Chuer DS eam ecauer If continued 
fraction expansion 
21-40 hos the second quotient eee 
41-60 he» the third quotient ZO 
61-80 ye ee OU rane ae keri O20 13 
N 1-20 Ayve7? the (4N-7)th quotient BZOeel3s 
Se iste chem coi 
21-40 Aung? 7 Oeee 4 P20 
41-60 Ayes? eN=o< 2M p20.13 
61-80 Nuys ¥N-4<2M D20.13 


Four quotients per card until 
2*M quotients have been input, 
where M is the system order. 
Assume this is the Lth card. 
The (L+1)th data card begins 
the second data set. 


Ze 





L+1 1-3 M=system order 3 


Ym 6 i inmor Cauerputleand K=2 for tS 
Cauer II inversion 


— 


The computer program has been written to handle up to 
20th order systems (M<20). If it is required to work with 
higher order systems, only one card change must be made. 


The specification statement is modified to read: 
REAL*8 A(N), BC(N), C(N)/N*0.Z2, D(N)/N*0./, DZERO 


where N is an integer no larger than 999. This restriction 


can be lifted by changing statement 2 to read: 
2 FORMAT( 2IR) 


where R is the mantissa of log, ){N). The REAL*8 in the 
specification statement indicates that all following varia- 
bles and arrays are real valued and in double precision. 
Modification to either single or extended precision would 
require changes in all format statements. If this is 
desireable, the user should consult references [25] and [26].: 
Execution time has shown to be less than .18 seconds 


for systems of order 10 or less. 
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-a oun, DD + 

REAL*8 A(20) ,8(20) ,C(20)/20*0./,0(20)/20*0./,DZERO 
eee READ IN SYSTEW ORDER, & WHICH CAUER INVERSION ... 


READ(5,1) MK 
FORMAT(213) 


eee READ IN QUOTIENTS FROM CONT. FRAC. EXPANSION wee 


READ(5,2) (A(T) ,B01),1=1_M) 
FORMAT (4D20.13) 


eee DETERMINE IF INVERSION IS CAUER I OR CAUER II wee 
IF(KeEQe1) GO TI 10 
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eoe INITIALIZATION oe. 
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APPENDIX C 


DETERMINATION OF WEIGHTING MATRIX 
(Q) FOR PRESCRIBED EIGENVALUES 


This FORTRAN IV program was used exclusively on the 
Naval Postgraduate School's IBM 360/67 digital computer 
and includes the associated job control language statements. 
The program consists of a main program and nine subroutine 


subprograms. The purpose of each subroutine is delineated 


below. 

SUBROUTINE Pesce rE. LON 

READ read in coefficients of numerator and denomi- 
nator polynomials of both transfer functions, 
and places each system in phase variable form. 

RAMAT determines the Routh array matrix and the 
transformation matrix, P. 

MULT PH multiples two matrices, Y and Z, and gives 
the resulting matrix YZ. 

HMAT RX determines the quotients of continued fraction 
expansion, H1(H2), and the state matrix in 
Gaver sl) tonm, HAIGH). 

POLYNM determines the product det (SI-A)xdet (SI+A" ) 

Mia LP computes the matrices G, T, and W as given 
in Chapter IV. ~ 

QTILDA computes the matrices, oF and Q from results 
of subroutine HELP. 

Send determines the_off diagonal elements, ree 
of the matrix Q. J 

Re be writes all two-dimensional matrices. 
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The input required has been reduced to a minimum. Multiple 


data sets are possible; and input as indicated below: 


Card Columns Description Format 
1. 1-3 N; the system order in 
2 1-16, ue Cseeqenelh. 2 atc tl | ¢ F16.8 
7-32, tHe denominator coefficients of 
33-48 the known transfer function, 


— for an Nth order system, L 
cards are required, where L 
=K+l and K is the integral part 


Chal oe 
L+2 1-16, PE Oa eZ ere tlt 5. the F16.8 
17-32, numérator coefficients of the 
33-48 known transfer function, L 
oe: cards required. 
Zine 2 1-16, eso waeee lis 2 ee Tl F16.8 
17-32, the denominator coefficients 
33-48, of the transfer function with 
Pa prescribed eigenvalues. 
olor 2 1-16, Por taltscaeennls the eG 
17-32, numerator coefficients of the 


33-48, transfer function with pre- 
rs scribed eigenvalues. 


for multiple data sets, repeat 
the same prodecure. Each data 
set requires 4Lt+l cards. 


This program has been written to accept systems up through 
20th order. To increase the capability of the program, only 
the dimension statements and the second continuation card of 
the equivalence statement require modification. The system 
order capability can be increased to 50. Beyond 50th order, 
the program requires an excessive amount of storage space 
(>510K bytes). Even this limitation is easily overcome by 


removing the four cards between statements 1000 and 1001. 
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Any other BEER one (i.e. , Single or extended pre- 
cision) require extensive changes to all subprograms. In 
GmemeGase, the user should consult [25] and/or [26]. It is 
recommended that an object deck or disk storage be used when 
available as compilation time is approximately 70-80% of 


total CPU time. 
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